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Abstract. Current distributions in thin film wires bent at different angles were investigated with 
a conformal mapping method. The technique of angles rounding using three parameters was sug- 
gested. This technique enabled to obtain a smooth line having a similar form with an arc of a circle 
and to avoid infinite current density in the angle. The dependency between current density and 
the radius of rounding was examined. 
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Introduction 

Current distributions in flat conductors are problems of interest because of printed circuit boards 
and on-chip devices development |14[ |4l llOj . Commonly these problems are solved numerically 
or using near- field measurements [5J [T]. In the event if there are various angles in conductors the 
numerical solutions have bad convergence near the angles because formally the current density in the 
angle is infinite. In real conductors there are not such problems because there are not absolutely sharp 
angles and a small rounding always exists. It is difficult to take this rounding into account numerically 
therefore analytical methods are required. Such methods were investigated by P. M. Hall O O [7] 
and L. N. Trefethen [13] but only the analytical solution for angle of 90° with rounding was found 
[5]. In this paper we have considered current distributions in thin film wires bent at different angles 
with rounding and suggested the technique of estimating an optimal radius of corner rounding in such 
wires. 



1. The conformal mapping for a conductor bent at arbitrary angle 

The geometry of a considered conductor is depicted in figure [T] a is an arbitrary angle. Only the case 
of h > k have been investigated (the other may be examined in a similar way). We have considered a 
direct current. This fact has enabled us to reduce Maxwell's equations to the Laplace equation for a 
scalar potential 

Aifi = 0. (1) 

This equation was solved for a complex potential W{z) = U{z) + iV{z) using the conformal 
mappings method; V{z) is a scalar potential and U{z) is a stream function. These two functions are 
bound with Cauchy-Riemann conditions therefore a boundary conditions can be written only for one 
of them. The boundary conditions requires the absence of a current flow through lateral boundaries 
of the conductor so the resulting boundary value problem is the following [3] 
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Ui^^y)\ix,y)en^ =Ui= const 
. Uix,y)\(^^^y)en2 ^U2 = const 

rii and Q2 are the upper and lower boundaries of the conductor respectively. 

To solve this problem we have obtained first the solution of Laplace equation in the upper 
complex half-plain with the boundary conditions corresponding a point charge at the origin. Than we 
have mapped this solution onto the considered domain. 

The required mapping corresponds the Schwarz-Cristoffel transformation of the following type 

/3 = a/iT. Undefined constants C and a have been determined from the conditions of boundary 
conformity and having the following view 

C ~ — (sin a — i cos a) 

TT 

where h, k and a are depicted in figure [T] 

The expression Q is representable in terms of elementary functions if 1 — /3 = P/Q where 
< P < Q; P and Q are integers [8]. From a physical point of view this condition docs not impose 
any restrictions on an angle because any irrational number can be approximated by the rational one 
with arbitrary high accuracy. Thus, ^ can be reduced to 

z = -{l + a)QcJ ,P(,g +a)(tQ-l) ^^ (3) 



where 



t = 



zi + a\Q 
zi - 



We have integrated the expression ^ for the angles of 60°, 30°, 120° and obtained the following 
results: 
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Figure 2. Distribution of a current density for a = 60° and h/k = 2. 



1. a ==60°: 



: C I ^ ln(i + ft) - ^ Ht" -bt + b^) - ln(i - 1) + ^ Ht" + i + 1) 



iv^ / -2i + 6 - i\/36\ is/l ^^f-2t-l-iy/i 



(4) 



where 



2. a = 30° 



h + k 



t = 



zi 



.V3 

-&3 



{7k - h) 



zi-1 
k 



-t + ibj ' 265 
/ 2t-l + 



V3bt 
VSbt 



1^ 



' 2t + l + iV^ \ 
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Figure 3. Distribution of a current density for a = 30° and h/k 

where 

C=A(i_,V3) 



= fc - /i^ - i Q + fc^/3^ 



1 



t = 



b- 



1 

h\ 5 
k 



a = 120° 



= C I ln{t + b) + ^ Ht^ -bt + 6^) - ln{t - 1) + ^ H*^ +t + l) 

2b ^ y-2t + b-iV3b ^ \^-2t-l- 

where 

C, = ^+^^i7k + h) 



t = 



zi 

zi-1 
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Figure 4. Distribution of a current density for a — 120° and h/k = 1.5. 

The lines of current calculated for these angles are shown in figures [2] ~ [4j 

The current strength J applied to the conductor supposed to be specified. The current distribu- 
tion for y — > — oo (fig. [2]|4]) is uniform therefore one can assume that current density at an infinite 
distance is j-oo = J/k. On the other hand, current density is proportional to the amplitude of complex 
potential W(zi) derivative |T^ : 



J{z) - A' 



dW{zi{z)) 



dz 



The complex potential of a point charge placed at the origin W{zi) ^ In(zi) [T^ therefore one can 
obtain the expression for the current density in terms of the implicit function zi{z) 

1-13 



J{z) = ^ 



dW{zi) dzi 



dzi dz 



= A 



1 f Zl{z) + I 



C \zi{z) - 1 



where A' and A are proportionality factors. 

y — > — oo and x = const correspond to value zi — rie^ 
values into ^ gives 



(7) 



The substitution of these 



= A 



whence it follows that 



A^ 



■ib 



J 



^A- 



(8) 



Finally, the distribution of current density is the following 

1-/3 



J 



1 / Zl{z)+i 



C \zi{z) - 1 



(9) 



The dependency zi(z) is given by the expression ([3| (or Q - Q) in implicit form and can be 
calculated numerically. The dependency between current density and the distance from the point B 
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Figure 5 . The dependency between scaled current density and scaled distance from 
the point B along the hne BE (fig. [T]) for the angles of 30°, 60° and 120°. / is the 
length of BE. 
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Figure 6. The points conformity for a conductor with the rounded angle. 



along the line BE (fig.[T]) is depicted in figure [5] The scaling is used in order to compare distributions 
for different angles. 

The point E of the complex plane z corresponds zi = 1 + i • (fig. [T]) . With relation to (|9| one 
can see that the current density in the angle is infinite. This result is nonphysical because in real wires 
there aren't absolutely sharp angles. 



2. Current distributions in wires with rounded angles 

To take the rounding into account we have used a technique similar to the one considered in [9] but 
we have used three new parameters instead of two. It was done in order to obtain smooths curve that 
looks like an arc of a circle. 
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We have modified the expression ^ in the following way 



zr = C 



1 fzi-l 



z\ \ zi + a 



1-/9 



dz\ 



^_fzi-l + 5i 
z\ 



dz\ + C7 



\_ / zi-l-(52 

z\\ zi + 5 



1-/3 



dzi (10) 



zi + a 

where constants (52 and 7 define the radius of rounding (fig. [6]). 

Constants a and C were defined in a similar way to the non-rounded case from the conditions of 
boundary conformity 

(sin a — icosa) 



C 



7r(l + 27) 



h\ 1-/3 r 1 + 7(1 - bxf-f^ + 7(1 - 82f-^ \ ^^f^ 



1 + 27 

To calculate the expression (10) we have supposed again that 1 — (3 = P/Q and denoted 

Zfl = /(l)+7/(l-<5i)+7/(l + ,52) 



where 



Changing variables 



I{m) = C 



1 / zi — m 
zi \ zi+a 



zi + a \ Q 



dzi. 



zi — m 



we have led /(m) to the rational integral 

I{m) = -{l + bl)QC 





fPffQ :hQ](fQ _ 1 



dt. 



(11) 



The obtained expression is formally equivalent to ([s]) therefore previously calculated expressions 
Q - ([6]) have been used to find the result. 

To determine the dependency between the radius of corner rounding p and parameters ^i, 82 
and 7 we have assumed that the rounding curve is an arc of a circle (fig. [?]). This assumption is 
correct for p ^ h,k. The placement of the point E corresponds to z(l) and the placement of the point 
E' corresponds to zj^(l). The lengths Ax and Ay were expressed through the radius p using simple 
geometric considerations 

(12) 
(13) 



Ax = Re[zfl(l)] - = p ( 1 - sin - 



Ay = yE - Im[z_R,(l)] = p ( 1 - sin 



cot 



The equivalence is inaccurate because the real curve isn't an arc of a circle. The coordinates of the 
point E are the following 



Xe — k 
h + k cos a 



(14) 
(15) 
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Figure 7. Schematic diagram of a rounded angle. 



The formulas (12 1 - (15) give two equations for three parameters. To obtain the third equation 

EH (fig. |6| . Finally, the system of equation for parameters Si , 



we have imposed the condition GE 
62 and 7 has been obtained 



Re[z^(l)] = k + p (l — sin 



lm[zR{l)] = - 



h + k cos a 



sma 



— p ( 1 — sin ■ 



cot 



(16) 



S2W 



\zr{1) - ZRil - Si)\^ = |zfl(l)-zj^(l 

The expression for current density for such mapping is defined in a similar way to ([9| and has 
the following view 



C [zi - ly-P + 7[0i - 1 + + 7[zi -1-62] 



1-/3 



where zi = zi{zr) 



The system (16) was solved numerically for angles of 120°, 60° and 30°. It was determined that 



the solution of it exists only for the point E' placed sufficiently close to the point E. That imposes 
the restriction on the value of p for every given angle. For instance, acceptable radii for the angle of 
120° should be less than 0.8fc for the considered case. For the angle of 30 ° the maximal radius is just 
0.012fc. Probably, bigger radii should be obtained using some other techniques. More detail approach 
to this problem is the object of further investigation. 

The dependencies between the density of current at the point E' normalised to j^o^ ^ and the 
radius p of angle's rounding normalized to the value of k are shown in figures |8) [3l and [12 



The current distributions in the conductors with the rounded angles are presented in figures |9j 
[IT] and [H 
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Figure 9. Distribution of current in the conductor with the rounded angle of 60° 
and h/k = 2. The rounding radius p = 0.02fc. 
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Figure 1 1 . Distribution of current in the conductor with the rounded angle of 30° 
and h/k = 2. The rounding radius p = O.Ollfc. 




Figure 13. Distribution of current in the conductor with the rounded angle of 120° 
and h/k = 1.5. The rounding radius p = O.Ik. 
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Conclusion 

The obtained results show that the consideration of current distribution in a thin film wires with 

absolutely sharp angles gives a nonphysical results, namely, an infinite current density in the angle. 
Therefore the technique of angle's rounding was suggested. Using of three parameters defining the 
rounding radius instead of two enables us to obtain a smooth line having a similar form with an arc 
of a circle. The obtained dependency between the extreme current density in a wire and a rounding 
radius can be used to estimate a radius of an angle's rounding required to prevent a destruction of a 
wire with an applied current. 
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